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ABSTRACT 

Models of planet formation and evolution predict that giant planets form efficiently in protoplanetary 
disks, that most of these migrate rapidly to the disk’s inner edge, and that, if the arriving planet’s 
mass is < Jupiter’s mass, it could remain stranded near that radius. We argue that such planets 
would be ingested by tidal interaction with the host star on a timescale < 1 Gyr, and that, in the 
case of a solar-type host, this would cause the stellar spin to approach the direction of the ingested 
planet’s orbital axis even if the two were initially highly misaligned. Primordially misaligned stars 
whose effective temperatures are > 6250 K cannot be realigned in this way because, in contrast with 
solar-type hosts, their angular momenta are typically higher than the orbital angular momentum of 
the ingested planet as a result of inefficient magnetic braking and of a comparatively large moment of 
inertia. Hot Jupiters located farther out from the star can contribute to this process, but their effect 
is weaker because the tidal interaction strength decreases rapidly with increasing semimajor axis. We 
demonstrate that, if ~ 50% of planetary systems harbored a stranded hot Jupiter, this scenario can 
in principle account for (1) the good alignment exhibited by planets around cool stars irrespective of 
the planet’s mass or orbital period, (2) the prevalence of misaligned planets around hot stars, (3) the 
apparent upper bound on the mass of hot Jupiters on retrograde orbits, and (4) the inverse correlation 
between stellar spin periods and hot-Jupiter masses. 

Subject headings: planets and satellites: dynamical evolution and stability — planet-star interactions 
protoplanetary disks — stars: rotation 


1. INTRODUCTION 

One of the puzzling findings in the study of exoplan¬ 
ets has been the difference in the properties of planets 
orbiting cool and hot stars. In particular, using mea¬ 
surements of the Rossiter-McLaughlin effect, it has been 
inferred that the apparent obliquity A (the angle be¬ 
tween the stellar spin and planet’s orbital angular mo¬ 
mentum vectors as projected on the sky) is typically 
small (A < 20°) for cool stars, but that stars of effec¬ 
tive temperature T e g > 6250 K exhibit a broad range of 
appa rent obliquities, reaching all the way up to A ~ 180° 
(e.g.. lWinn et al.lfeOlOHAlbrecht et al.ll2012f) . It was fur¬ 
ther found that the masses M p of planets on apparent 
retrograde orb its (A > 90°) are < 3Mj (where Mj is 
Jupiter’s inass: lH6brard et aLll201lH . and that the stellar 
rotation periods P* of cool stars t hat host a hot Jupiter 
(HJ) decrease with increasing M p dDawso il l20ll . 

Several models have been advanced to explain the 
obliquity dichotomy, but so far none can account for all 
the relevant observations. Interpretations based solely 
on di fferences in the stellar properties (e.g.- lRogers et all 
12013) cannot address the inferred dependence of the 
obliquity on M p and the detection of highly misaligned 
planets around cool stars. In view of the fact that T e ff ~ 
6250 K corresponds to the temperature above which the 
size of the outer convective zone of a main-sequence (MS) 
F star shrinks rapidly, it was suggested that close-in plan¬ 
ets in both cool and hot stars are initially distributed 
over the entire angular range (0 — 180 degrees), but that 
only in cool stars where a substantial convective enve¬ 
lope is present can a sufficiently massive close-in planet 
(M p > lMj) realign the star through tidal interaction 
(e.g.. IWinn et all 120101 : lAlbrecht et al.ll2012D . The host 


stars are also subject to magnetic braking, which declines 
strongly above the same transition temperature (corre¬ 
sponding to the br eak in the Kr aft curve; IKra 1lM3), 
and it was argued (|Dawsonll2014f) that this, rather than 
the tidal dissipation efficiency, is the main factor under¬ 
lying the difference in obliquity properties between cool 
and hot stars. Both variants of this scenario, however, 
face the conundrum that an HJ undergoing equilibrium 
tidal interaction with its host star would spiral in and be 
ingested on a timescale that is comparable to the realign¬ 
ment time. A possible way out of this difficulty is to ap¬ 
peal to the {10} component of the dynamical tide, which 
could in principle significantly reduce the alignm ent time 
without affecting the ingestion time (Lai 12012} FI How¬ 
ever, even though t his model can be used to account for 
individual systems (IValsecchi fc Rasioll2014f) . it remains 
unclear whether it can explain the overall A distribution 
of HJs an d the manifested d i fference between cool and 
hot stars dRogers fc Linll2013i : IXue et al.ll2014l l. Further¬ 
more, even the basic tidal interaction interpretation of 
the obliquity dichoto my has now been c alled into ques¬ 
tion by the results of lMazeh et al.1 (12015I1 . who analyzed 
the rotational photometric modulations of a large sam¬ 
ple of Kepler sources and inferred that (1) the conclusion 
that planets around cool stars are well aligned, and those 
around hot stars are not, is general and not restricted just 
to HJs; and (2) the low obliquity of planets around cool 
stars extends to orbital periods P OI b that are a factor of 
~ 10 larger than the maximum value (< 4 days) for ro¬ 
bust tidal interaction between an HJ and a > 1 Gyr-old 

1 An alternativ e possibility, that only the out er convective la yer 
llWinn et aUl201 t)j — or even just a part of it (IDawso lESD - 
partakes in the realignment process, is hard to justify on either 
theoretical or observational grounds. 
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G or F star. 

In this Letter we propose to address these apparent 
difficulties and account for many of the observed differ¬ 
ences between the properties of planets in cool and hot 
stars by postulating that, in addition to the tidal inter¬ 
action with existing close-in planets, a large fraction of 
the stellar hosts — both cool and hot — ingest a hot 
Jupiter early on in their evolution. This proposal is mo¬ 
tivated by the ex pectatio n that a la rge fraction (up to 
80% according to lTrilling et al.lf2002 ) of solar-type stars 
possess giant planets during their pre-MS phase, and that 
a large fraction of the giant planets that form in a pro¬ 
toplanetary disk on scales < 5 AU migrate clo se to their 
host star before the disk is dispersed (e.g., Ilda fc Lffil 
I2004T) . Numerical simulations incorporatin g an N-body 
code and a ID a-viscosity disk model dThommes et al.l 
I2008T) demonstrated that this behavior can be expected 
for disks with a < 0.01 that are sufficiently massive. The 
inward planet migration is likely stopped by the strong 
(> 1 kG) protostellar magnetic field that tru ncates the 
disk a t a radius (r; n ) of a few stellar radii fe.g.. lLin et al.l 
fl996h . Gravitational interaction with the disk causes a 
planet reaching n n to penetrate into the magnetospheric 
cavity and, if it is massive enough, to undergo eccentric¬ 
ity excitation that can rapid ly lead to a collision with 
the star (e.g.. lRice et alJl2008lh It was, however, inferred 
that if M p is sufficiently small (< 1 Mj), the planet would 
remain stranded at a distance where its orbital period is 
~ 0.5 of that at r; n until well after the gas disk disap¬ 
pears (on a timescale of ~ 10 6 — 10 7 yr). In our proposed 
scenario, the primordial disk orientations span a broad 
angular range that is reflected in the orbital orientations 
of the stranded planets. When the latter are ingested 
by tidal interaction with the host star (on a timescale 
< 1 Gyr), the absorbed angular momentum is sufficient 
to align a solar-mass star in that general direction, but 
not an MS star with T e s > 6250 K. This is because cool 
stars have significantly lower angular momenta at the 
time of ingestion than hot stars as a result of a more effi¬ 
cient magnetic braking process and of a lower moment of 
inertia. Given the proximity of the stranded HJs (SHJs) 
to their host stars (P or b,SHj < 2 days), they can be ex¬ 
pected to have been ingested by the time their parent 
planetary systems are observed; however, giant planets 
farther out can continue to interact with their host stars 
and potentially affect their measured obliquities. In our 
simplified formulation, we model the SHJs using as pa¬ 
rameters their characteristic mass Mshj and the fraction 
p of systems that initially harbored an SHJ. By compar¬ 
ing the predictions of this model with the observational 
data, we infer Mshj < 1 Mj and p ~ 0.5. 

The potential effect of tidally induced ingestion of HJs 
on the observed properties of plane tary systems ha s been 
recognized before. In particular, I Jackson et al.1 (120091 1 
suggested that this process could account for the ob- 
served orbital di s tribut ion of close-in planets, whereas 
iTeitler fc Konigll (12014H proposed that the spinup in¬ 
duced by the deposition of a swallowed planet’s orbital 
angular momentum in the host’s envelope can explain the 
observed dearth of close-in planets around fast-rotating 
stars. The current proposal further extends this scenario 
by incorporating not only the distribution of observed 
planets but also a putative population of HJs that were 
ingested after being stranded near the inner edges of their 


associated protoplanetary disks. 

2. MODELING APPROACH 

Using Monte Carlo simulations, we follow the temporal 
evolution of the stellar and orbital angular momenta of 
10 6 systems in the context of an equilibrium-tide model. 
The Cartesian coordinate frame [x, y, z) is chosen so 
that the orbital angular momentum (L) always points 
along the z axis and the stellar one (S) always lies in 
the y-z plane. We assume circular orbits and neglect the 
precession of L and S around the total angular momen¬ 
tum vector 0 The set of equations employed to model 
the tidal interaction between a star (subscript *) and a 
planet (subscript p) is: 

dL L 12*3, „ L 12*3, „ L i 12*- \ „ , 

dt 2 t ,] f2 or b 2r t ] !2 or b Ti y !2 0 rb J 

dS L 12*3, - dj 12*3, „ , dL i I2* z \ „ 

— =- x - —y H-I 1 —- 1 z — T . 

dt 27b I2 0 rb 2 t,| I2 0 rb Tl \ Uorb / 

( 2 ) 

where f2 or b = f2 or b z and $2* are the orbital and stel¬ 
lar angular velocities, respectively (related to the orbital 
and stellar rotation periods through P or b = 27r/f2 or b and 
P* = 27r/f2*), T = TS is the magnetic braking torque, 
and 7d is the nominal tidal damping time, which, neglect¬ 
ing dissipation inside the planet, is given by 

4Q'*/a\ 5 M* 1 

Td 9 \Rj M p f2 orb ’ [) 

In Equation ([3]), Q* is the tidal quality factor (taken 
to be a constant, for simplicity), a is the semimajor 
axis, whereas R and M denote an object’s radius and 
mass, respectively. We model the star as a uniformly 
rotating body (S = 7*f2*) with a moment of inertia 
/* « 0.0 6 M*P^, and adop t for T the expression pre¬ 
sented bv lMatt et al„ (2015t h which covers both the sat¬ 
urated and the unsaturated regimes and incorporates as a 
key variable the convective turn over timescale as a func - 
tion of T c s (3300 < P e ff < 7000: iCranmer fe Saadl20lTlh 
The nominal magnetic braking time is r m b = |SyTj. 

The equations are solved in a two-step process: we first 
integrate for a small At employing a 4th-order Runge- 
Kutta method and then rotate the vectors so that the 
updated L is along z and the updated S lies in the y-z 
plane. When considering a multiplanet system, we solve 
the equations for the different planets in sequence at each 
time step, neglecting possible interactions between them 
and assuming that their orbits remain coplanar. 

We consider two sets of models based on the type 
of host star: G or F; see Table [1] for a listing of the 
adopted parameters. For each stellar type, we evolve 
5 x 10 5 systems, in each case considering simultaneously 
the SHJs (parametrized by Mg H j and p\ see Section Q]) 
as well as 5 additional planets 0 for which we randomly 
pick the orbital periods, planet radii, a nd system age ac- 
cording to the procedure described in ITeitler fe Konigll 

2 The associated torques act along the x axis and thus do not 
affect either the orbital separation or the alignment. 

3 Our results are not sensitive to the choice of the number of 
additional planets because large, close-in planets — which, accord¬ 
ing to Equation ©. exert the strongest tidal effect — are a rare 
occurrence. 
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TABLE 1 

Model Parameters 


Parameter 

Value (G / F) 

Distribution 

T eS [K] 

5500 / 6400 


M * [Mq] 

1.0 / 1.3 


R* [R®] 

1.0 / 1.3 


R* ini [days] 

5-10 / 2.5-5 a 

Uniform 

Age [Gyr] 

1-8 / l-4 b 

Fit to data 

Number ot planets 

5 


Rorb [days] 

0.5-50° 

/(In P) oc P 0 ' 54 

R p [R®] for P olh < 7 d 

2-20° 

/(lnRp)ocRp 109 

R p [R®] for P orb > 7 d 

2-20° 

/(InRp) OCR" 2 ' 31 

M p [Mq] for R p < 9R® 

4-8 l d 

(R p / Rq) 2 Mq 

M p [M®] for R p > 9R® 

56-5620 d 

Fit to data 

V’ini [degrees] (random) 

0-180 

/(b) oc sin(b) 

p (SHJ fraction) 

0.5 


-^SHJ [M, j] 

0.6 


Rshj [Rj] 

1.0 


Rorb,SHj [days] 

2.0 



References. — a iMeihom ct at (2011); b IWalkowicz fc Basnl 
l20T3) ; c lYondinl OoTilh d I Weiss et all <20l3) . 


Note. — The subscript “ini” indicates an initial value. 

(120141) PI As in that work, we adopt the approxima¬ 
tion that the entire planetary orbital angular momen¬ 
tum is added to S when the planet reaches the Roche 
limit or « 1.5 (3 M*/Mr,) 1 ' 3 R p, and we neglect st ellar 
evolution effects (but se e, e.g., IValsecchi et al.l[20l5l and 
lYalsecchi fe Rasiol 120141 ). The effective value of ()[ and 
its ph ysical basis are still open questions (e.g., lOgilvid 
I2014T) . Note, however, that the apparent break at P or b > 
3 days in the peri od distr i butio n of planets around solar- 
mass stars (e.g., lYoudinl [20111) . if due to tidal interac¬ 
tion is consistent with Q' t rs 10 4 5 6 (e.g.. iTeitler fc Konigll 
[Ml. Using the data from exoplanets.org, we find 
that F-star systems manifest a similar break at the same 
location. We therefore adopt Q ' ~ 10 6 for both our cool 
and hot fiducial stars. 

The initial spin-orbit angle ■i/'ini is taken to be random, 
f(ijj ini) oc sin(^i n i), which correspond s to a flat distribu¬ 
tion for the projected angle A ini (see iFabrvckv fc Winnl 
12009T ). Our choice of P 0 rb,SHj ~ 2 days is motivated 
by th e characteristic period (~ 7 days; iGallet &; Bouvierl 
I2013D of disk-accreting protostars, which can be inter¬ 
preted i n ter ms of the magnetic disk-locking model (e.g., 
iKonigll ri991 1 In this picture, the disk truncation ra¬ 
dius is located at r- m ~ 0.7 r co (with r co defined by 
Uorb(T’co) = U*) if the magnetic in teractio n is dorn i- 
nated by the dipolar field component dLong et al.ll2005D . 
The choice of P or b,SHj then follows from the results of 
iRice et all (120081) . who inferred that it is ~ 0.5 P O rb(Un)0 

3. RESULTS 

The results for the final configurations of the modeled 
planetary systems are presented in identical formats in 

4 In relati ng planet masses to pla net radii, we generalized the 
treatment in ITeitler &; Konigll ([20141 ) by employing a power-law 
relationship only for R p < 9 R& and using an empirical fit to the 
data for R p > 9 R@; see Table 111 

5 Using the numerical results of ILong et al.l (120051) and 
the characteristic radius (~ 2 R@) and mass accretion rate (~ 

10 —8 Mq yr —x ) for classical T Tauri stars, we infer a stellar dipo¬ 
lar field strength of ~ 2.5 kG, which is consistent with observations 


TABLE 2 
SHJ Ingestion 


Parameter 

Init. 

G star 
Roche 

Ingest. 

Init. 

F star 
Roche 

Ingest. 

t [Gyr] 

0.0 

0.69 

0.69 

0.0 

0.35 

0.35 

a [AU] 

0.03 

0.01 

— 

0.03 

0.01 

— 

Rorb [days] 

2 

0.5 

— 

2 

0.5 

— 

R» [days] 

7.5 

12.9 

5.0 

3.8 

4.5 

4.2 

Anb [Gyr] 

0.21 

0.61 

0.09 

1.96 

2.83 

2.46 

|Rshj/S| 

1.59 

1.71 

— 

0.43 

0.33 

— 

-0 [degrees] 

100 

35 

13 

100 

86 

69 


Figures Q] and [2] for our representative G and F stars, 
respectively. The top three panels show the probabil¬ 
ity distribution function of the projected spin-orbit an¬ 
gle for the entire population (left) and separately for 
its short-period (P or b < 5 days, middle) and long-period 
(Porb > 5 days, right) components. These panels demon¬ 
strate that our model can ac c ount for the basic trends un¬ 
covered by the lMazeh et al.l (120151) study — that planets 
around cool stars are well aligned with their host’s spin 
irrespective of the planet’s size or orbital period, but that 
planets around hot stars are not. To q uantif y the degree 
of alignment, we follow iMazeh et al.l (120151) and evalu¬ 
ate the ratio q = (sin(i*)) /(7 t/ 4), where i* G [0, 7r/2] 
is the angle between the stellar rotation axis and the 
line of sight, the angle brackets denote an average over 
the distribution, and the normalization is by the value 
of (sin(i*)) for a random orientation of S on the plane 
of the sky. This ratio is therefore 1 when there is no 
preferred spin-orbit orientation (as in our adopted ini¬ 
tial condition), whereas it is q = 1 /(tt/ 4) ~ 1.273 when 
S is normal to the line of sight (which, in the case of a 
transiting planet, likely corresponds to A « 0). For our 
G-host model we obtain q = 1.145, which corresponds 
to having the ~ 50% of systems that ingested an SHJ 
attain good alignment even as the rest remain close to 
their initial random orientations. For the F-host systems 
we find q = 1.030, indicating that in this case the initial 
random distribu tion is o nly sl ightly modified. We note 
from Figure 1 of IMazeh et all (120151 ) that, while the me¬ 
dian value of q for stars with T e g < 5000 K is consistent 
with perfect alignment, this ratio decreases toward 1 in 
the G-star range (T e g ~ 5300 — 6000 K) that we have 
modeled^ 

The bottom three panels of Figures []] and [2] confront 
our model predictions for HJs (defined here by M p > 
0.5 Mj and P or b < 10 days) with the observational data 
(as of June 2015) obtained from exoplanets . org, which 
consist of Nq = 36 entries for cool stars (T e g < 6250 K) 
and Np = 19 entries for hot ones (6250 K < T e g < 
7000 K). The left panel shows the planet counts as a 
function of A, with the model predictions obtained by 
selecting 1000 samples from our simulations, each con¬ 
sisting of Nq (Figure [TJ or Np (Figure [2j) systems, and 
plotting the average counts per bin (bars) as well as the 
range of ±1 standard deviation (shaded areas). It is seen 
that the model is in excellent agreement with the data, 

6 IMazeh et alj (120151 1 inferred q < 1 for their hot-star sample, a 
value that cannot be explained in the context of our basic model. 
Note in this connection that, as discussed by these authors, the 
observational challenges that need to be overcome in deriving the 
hot-star result are more severe than those encountered in the cool- 
star case. 
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with the only > la deviation exhibited by the bin en¬ 
compassing A = 180° Q 

The key feature that enables our model to account for 
the aforementioned observations is the ingestion, over 
time intervals (fingest) that are shorter than the ages of 
the observed systems, of SHJs with values of |Tshj/S'| 
that are > 1 for G stars and < 1 for F stars (see Ta¬ 
ble [2]). The low obliquity attained by initially misaligned 
G stars at ti ng est is due mainly to the efficient reduc¬ 
tion of S through magnetic braking over the time in¬ 
terval /-ingest.5 by Contrast, T m b(/ £ /-ingest.) A /ingest Icr 
F stars. The strongest constraints on the parameters 
of the postulated SHJ population are provided by the 
observed A distributions of HJs for the G- and F-star 
systems. In particular, by comparing the number counts 
of well-aligned HJs in solar-type stars (corresponding to 
the first bin in the lower left panel of FigurcQ]) with those 
in the A > 20° bins, one infers p ~ 0.5: the value of this 
parameter can be expected to reflect both the forma¬ 
tion rate and radial transport properties of HJs in pro¬ 
toplanetary disks and the fraction of systems in which 
disk truncation is eith er short-lived or inefficient (e.g., 
iHerbst fe Mundtl 1200511 . The corresponding data for F 
stars (lower left panel of Figurc[2j in turn imply an upper 
bound (~ lMj) on Mshj: If the SHJ mass were mea¬ 
surably larger, there would be significantly more aligned 
HJs detected in these systems. Interestingly, this upper 
bound is consistent with the li mit obtained i n the SHJ 
formation scenario studied bv lRice et al.1 (I2008fl . We em¬ 
phasize that the qualitative results of our model are not 
sensitive to the exact values of these (or any of the other) 
parameters. Additional data would, however, be useful 
for better constraining these values. 

The host stars also interact tidally with HJs that are 
still orbiting (and thus remain observable), which typi¬ 
cally have P or b values of a few days. While the impact 
of this interaction is not as pronounced as that of inges¬ 
tion, it can still affect the observed properties of these 
systems. This is illustrated in the middle and right bot¬ 
tom panels of Figures CD and O which display, respec¬ 
tively, the dependence of A and of the projected stellar 
angular velocity on the HJ mass for a random selection 
of 100 simulated systems. It was already recognized be¬ 
fore that the tidal interaction scenario is consistent with 
the apparent dearth of massive planets with either ret¬ 
rograde orbits or high-P* hostsO This follows from the 
inverse dependence of rj (Equation (|3jl) on M p , which 
implies that more massive HJs should be more efficient 
at realigning and spinning up their hosts. The displayed 
results demonstrate that our model also broadly repro¬ 
duces these trends quantitatively, for both G and F stars. 


relies on the orbital planes of the ingested SHJ and of 
any remaining distant planet having roughly the same 
orientation, which is consistent with the picture of gi¬ 
ant planets forming in a nearly planar disk and mi¬ 
grating to the vicinity of the host star. However, this 
scenario is not consistent with interpretations of pri¬ 
mordial misalig nment i n terms of ei t her planet-planet 
scattering (e.g., iBeauge fe Nesvornvl 120121) or a Kozai- 
Lido v-type i nteraction with a misaligned companion 
(e.g., iFabrvckv fc Tremaine! 12007! ). as these mechanisms 
typically result in different orientations for the orbit 
of an HJ and those of more distant planets. There 
have already been several proposals in the literature for 
forming protoplanetary disks that are misaligned with 
the protostellar spin axis, including a warping torqu e 
exerted by the stellar magnetic field (lLai et al.l 12011 1: 
a gravitatio nal torque exerted by a misaligned com¬ 
panion dBatvginl 120121) , possibly amplified by a reso¬ 
nance between the disk-torquing frequ ency and the (disk- 
drive n) stellar precession frequency dBatvgin fe Adams! 
2013T): a combination of th e above two torques ( Laill20l4 


Spalding fc Batygin! 1 201 111 : and accretion from a turbu¬ 


lent i nterstellar medium ()Bate et al .1120 1 Ol iFielding et al.1 
120151 ). However, these models have not yet been 
fully evaluated in light of the observed A distrib ution 
of transiting planets (e.g., iCrida fc Batvginl 120141) and 
other observational constraints (e.g.. lWatson et alJl20Tlt 
iGreaves et al .1 fjoTjl) . In this work we adopted a com¬ 
pletely random distribution of obliquities for simplicity, 
but more data and improved modeling might be able to 
constrain in the context of this scenario. Future 

steps toward refining our model could include taking ac¬ 
count of stellar evolution and of the effect of dynamical 
tides, incorporating an inner boundary condition that 
represents the stellar magnetosphere into global mod¬ 
els of planet formation and evolution to better predict 
the properties of SHJs, and reconciling the SHJ forma¬ 
tion picture with the observed distribution of lower-mas s 
planets (e.g., iMandell et aDl2007t [Fogg &; Nelsonll2007l) . 

The clearest observational prediction of this model is 
the occurrence of highly misaligned HJs with P or b < 
2 days around very young stars, but finding such ob¬ 
jects could be rather challenging fe.g.. lMiller et al. 1120081 
iCrockett et al.l 120121) . There is at present only one 
claimed detection of a short-period (P or b ~ 0.45 days), 
Jupiter-ma ss planet around a pre -MS (~ 3Myr), low- 
mass star (Ivan Evken et al.l l2012f ). and it has, in fact, 
been inferred to be highly misaligned (A ~ 70°); how¬ 
ever, the true nat ure of this object is still being debated 
(e.g., iCiardi et al.ll2015l and references therein). 


4. DISCUSSION 

Our interpretation of the good ali gnment inferr ed for 
long-period planets around cool stars ( Mazeh et al.ll2015h 
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Fig. 1.— System properties in the cool-star model, as functions of either the projected obliquity A or the projected stellar angular 
velocity. The top panels include the entire range of planet masses and show the A probability distribution for all the planets (top left) and 
separately for those with P Q rb < 5days (top middle) and P Q rb > 5days (top right). The bottom panels show results for hot Jupiters only, 
along with observational data taken from exoplanets.org: the lower left displays planet number counts and the other two the planet mass. 



A [degrees] 


A [degrees] 


A [degrees] 





Fig. 2.— Same as Fig. [T] but for the hot-star model. 
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